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ABSTRACT 

We report the first Very Long Baseline Array (VLBA) observations of 43 GHz 
v—1, J =1-0 SiO masers in the circumstellar envelope of the M-type semi- regular 
supergiant variable star AH Sco at 2 epochs separated by 12 days in March 2004. 
These high-resolution VLBA images reveal that the distribution of SiO masers 
is roughly on a persistent elliptical ring with the lengths of the major and minor 
axes of about 18.5 and 15.8 mas, respectively, along a position angle of 150°. And 
the red-shifted masers are found to be slightly closer to the central star than the 
blue-shifted masers. The line-of-sight velocity structure of the SiO masers shows 
that with respect to the systemic velocity of —6.8 km s _1 the higher velocity 
features are closer to the star, which can be well explained by the simple outflow 
or infall without rotation kinematics of SiO masers around AH Sco. Study of 
proper motions of 59 matched features between two epochs clearly indicates that 
the SiO maser shell around AH Sco was undergoing an overall contraction to the 
star at a velocity of «13 km s _1 at a distance of 2.26 kpc to AH Sco. Our 3- 
dimensional maser kinematics model further suggests that such an inward motion 
is very likely due to the gravitation of the central star. The distance to AH Sco 
of 2.26±0.19 kpc obtained from the 3-dimensional kinematics model fitting is 
consistent with its kinematic distance of 2.0 kpc. 

Subject headings: circumstellar matter — masers — stars: individual (AH Sco) 
— stars: kinematics 



1. Introduction 

AH Scorpii (AH Sco) is a semi-regular variable with an optical period of 714 days 
(Kukarkin et al. 1969) and a spectral type of M5Ia-Iab (Humphreys 1974). The systemic 
velocity of AH Sco is estimated to be about —7 and —3 km s _1 based on the observations 
of OH maser (Baudry, Le Squeren & Lepine 1977) and H 2 maser (Lepine, Pase de Barros 
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& Gammon 1976), respectively. The distance to AH Sco remains uncertain. A photometric 
distance of 4.6 kpc has been derived by Humphreys & Ney (1974) based on their infrared 
data. However, Baudry, Le Squeren & Lepine (1977) led to a photometric distance of 2.6 
kpc under an assumption of an absolute visual magnitude of —5.8 for lab stars, whereas its 
kinematic distance was about 1.5 ~ 2.0 kpc for its systemic velocity —5.5 ~ —7.5 km s -1 . 

Late type stars often exhibit circumstellar maser emission in molecules OH, H 2 0, and 
SiO. The supergiant variable AH Sco is such a star that has been detected strong maser 
emission with single-dish in all three species (e.g. Lepine, Pase de Barros & Gammon 1976; 
Baudry, Le Squeren & Lepine 1977; Balister et al. 1977, Gomez Balboa & Lepine 1986). The 
interferometric observations of these masers would be useful in determining the structure and 
kinematics of the circumstellar envelop (CSE) and understanding the physical circumstance 
and mass loss procedure for this supergiant variable. Unfortunately, there has been so far 
no any published interferometric map of OH, H 2 and SiO masers toward this source. 

Especially, SiO masers provide a good probe of the morphology of CSE and kinematics 
of gas in the extended atmosphere which is a complex region located between the photosphere 
and the inner dust formation shell. Previous VLBI experiments have demonstrated ringlike 
configurations (e.g. Diamond et al. 1994; Greenhill et al. 1995; Boboltz, Diamond & 
Kemball 1997; Yi et al. 2005, Chen et al. 2006), or elliptical distributions (e.g. Boboltz & 
Marvel 2000; Sanchez et al. 2002; Boboltz & Diamond 2005) of SiO masers, and also revealed 
complex kinematics in SiO maser regions, e.g. contraction and expansion at the different 
phase of stellar pulsation (Boboltz, Diamond & Kemball 1997; Diamond & Kemball 2003; 
Chen et al. 2006) and even rotation (Boboltz & Marvel 2000; Hollis et al. 2001; Sanchez et al. 
2002; Cotton et al. 2004; Boboltz & Diamond 2005). The ringlike or elliptical distribution 
that is assumed to be centered at the stellar position with a radius of 2-4 R* suggests that 
SiO masers are amplified in tangential rather than radial path. 

In this paper, we present the first VLBI maps of SiO maser emission toward AH Sco 
observed at two epochs separated by 12 days in March 2004. The observations and data 
reduction are described in § 2; results and discussions are presented in § 3, followed by 
conclusions in § 4. 

2. Observations and data reduction 

The observations of the v =1, J=l-0 SiO maser emission toward AH Sco (a =17Hl m 16.98 s , 
5 = -32°19'31.2", J2000) were performed at two epochs on March 8, 2004 (hereafter epoch 
A) and March 20, 2004 (hereafter epoch B) using the 10 stations of the Very Long Baseline 



- 3- 



Array (VLBA) of the NRAC0. A reference frequency of 43.122027 GHz was adopted for the 
v—1, J =1-0 SiO transition. The data were recorded in left circular polarization in an 8 
MHz band and correlated with the FX correlator in Socorro, New Mexico. The correlator 
output data had 256 spectral channels, corresponding to a velocity resolution of 0.22 km 
s -1 . The system temperatures and sensitivities were on the order of 150 K and 11 Jy K -1 , 
respectively, for both epochs. 

For the data reduction, we followed the standard procedure for VLBA spectral line 
observations using the Astronomical Image Processing System (AIPS) package. The band- 
pass response was determined from scans on the continuum calibrator (NRAO530). The 
amplitude calibration was achieved using the total-power spectra of AH Sco based on the 
"template spectrum" method. The template spectra at each epoch were obtained from the 
Mauna Kea (MK) station at a high elevation. A zenith opacity of about 0.05 at MK station 
estimated from the variation of system temperature with zenith angle was applied to correct 
the atmosphere absorption for each epoch. Residual group delays determined from a fringe 
fitting to the continuum calibrator were applied to the spectral line data. Residual fringe- 
rates were obtained by performing a fringe-fitting on a reference channel (at Vlsr = —9.4 
km s -1 ), which has a relatively simple structure in the maser emission. An iterative self- 
calibration on the reference channel was performed to remove any structural phase. The 
solutions of fringe-fitting and self-calibration were then applied to the whole spectral line 
data. Image cubes were produced for all the velocity channels between 14 and —18 km 
s -1 with a synthesized beam of 0.69 mas x 0.20 mas at a position angle —6°. Off-source 
rms noise (cr rms ) in channel maps ranges from 20 mJy beam -1 in the maps with weak or 
no maser emission to 50 mJy beam -1 in the maps containing strong maser emission. The 
flux densities, and positions in right ascension (R.A.) and declination (Dec.) for all emission 
components with intensity above 8 cr rms in each channel maps were determined by fitting a 
two-dimension Gaussian brightness distribution using the AIPS task SAD. Errors in R.A. 
and Dec. obtained from above fitting procedure range from 1 /ias for components with high 
SNR, to 168 /ias for components with low SNR, and the typical uncertainty of the fitted 
position of maser components was smaller than 10 //as. 

The remaining analysis of the maser component identifications was performed outside 
of the AIPS package. As described in our previous work (Chen et al. 2006), a maser spot 
is a single velocity component of the maser emission in each velocity channel map; a maser 
feature is a group of the maser spots within a small region in both space and Doppler velocity, 
typically 1 AU and 1 km s -1 , and is expected to be a physical feature consisting of a single 
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gas clump. In order to study the characteristics of SiO masers, it is necessary to identify 
maser features for each epoch. The maser spots in different channels were deemed as the 
same feature according to the criterion that these spots appear in at least three adjacent 
channels and lie within an angular separation of 0.5 mas. Finally, 82 and 87 maser features 
were identified for epochs A and B, respectively. 

3. Results and Discussions 

3.1. The spatial structure of the SiO masers 

The full lists of parameters for each identified feature are given in Tables 1 and 2, for 
epochs A and B, respectively. We fit a Gaussian curve to the velocity profile of a feature 
containing at least four spots to determine Vlsr at the peak of velocity profile. For some 
features which can not be well represented by a Gaussian profile (labelled by a "*" in Tables 
1 and 2), intensity weighted mean Vlsr was adopted. The velocity range across the feature 
Am, is defined to be the difference between the maximal and minimal velocities of the spots 
in the feature. Feature positions (x, y) in R.A. and Dec. were determined from an intensity 
weighted average over maser spots in the feature. The uncertainty (a x ,a y ) of a feature 
position was defined as squared root of the square sum of (1) the mean spot distance from 
the defined feature position and (2) the mean measurement error of the spot positions. The 
weights proportional to the intensity of the spot were applied in the uncertainty estimation. 
The typical position uncertainties of features are 0.01 mas and 0.02 mas for R.A. and Dec, 
respectively. The positions are measured with respect to the reference feature at (0, 0) for 
aligning the maps in the two epochs (labelled by "R" in Tables 1 and 2; see Sect. 3.2.2). The 
distance of a maser feature, r, is measured with respect to the fitted center obtained from 
the ellipse model fitting to the distribution of maser features (see below). The flux density 
of the brightest spot in each feature was deemed as its peak flux density, P. 

In two top panels of Figure 1, we compare the total power imaged by the VLB A (open 
circle) to the total power (solid line) obtained from the MK antenna for each epoch. The 
total power imaged by VLBA is obtained by summing all fitted flux of spots belonging to 
features. The fractional power representing the ratio of the total flux imaged by VLBA 
to the total flux of maser emission is shown in two bottom panels. The fraction is mostly 
between 0.4 and 0.8. That is, on average about ~60% of the total luminosity of masers 
was detected in our observations. Actually, such a fractional power reflects the degree of 
extension of the maser emission. If the apparent sizes of maser components are larger than 
that of the interferometric beam, the maser would be partly resolved (i.e. the fractional 
power is less than 1) and the fractional power should decrease with the increase of apparent 
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size of masers. The typical size of maser spots estimated from geometric mean of sizes of the 
major and minor axes of the spots, which were obtained by fitting to an elliptical Gaussian 
brightness distribution in the CLEAN map, is 0.5 mas. This typical scale size is slightly 
larger than the geometric mean of the VLB A beam of 0.4 mas. Thus ~40% missing flux in 
our map is mainly due to the high spatial resolution of the interferometric array. 

Figure 2 shows the distributions of maser features toward AH Sco for the two epochs. 
These high resolution VLBI images reveal a persistent elliptical structure of SiO masers 
around AH Sco druing an interval of 12 days. We characterized this morphology by per- 
forming a least-squares fit of an ellipse to the distribution of masers weighted by the flux 
density of each feature for each of two epochs. The best-fitting ellipses and the ellipse centers 
are also shown in Figure 2. The lengths of the major and minor axes were found to be 18.6 
and 15.7 mas for epoch A, and 18.4 and 15.9 mas for epoch B, respectively, with the major 
axis of the ellipse oriented similarly at ~ 150° at both epochs. And the fitted centers of the 
elliptical distributions are almost the same at both epochs of —7.9 and 5.8 mas in R.A. and 
Dec, respectively. At a distance to AH Sco of 2.26 kpc (see Sect. 3.2.3), the distribution of 
SiO masers corresponds to about 42 x 35 AU for both two epochs. However, an axial ratio 
of 1.18 suggests that the distribution of maser features around AH Sco can be viewed ap- 
proximately as the ringlike structure with an average diameter of about 17.2 mas (obtained 
from the geometric average of the major and minor axes for both epochs). 

The simultaneous near-infrared interferometry and radio interferometry imaging of cir- 
cumstellar SiO maser for late type stars have been done recently (e.g. Boboltz & Wittkowski 
2005; Wittkowski et al. 2007; Cotton et al. 2004, 2006). These observations reveal that the 
ratio of the maser ring radius to the photospheric radius of the central star is about 1.5 — 4.0. 
Unfortunately, there has been no any published photospheric radius measurement for AH 
Sco. Thus we can not directly compare SiO maser ring radius with stellar radius for AH Sco. 

We also notice that the red-shifted SiO masers lie slightly closer to the center than the 
blue-shifted masers (see Figure 2). This can be seen more clearly in Figure 3, showing the 
maser feature distance from the fitted center versus its line-of-sight (LOS) velocity (see Sect. 
3.2.1). This phenomenon seems to be explained under the assumption that maser gas was 
undergoing infall to the central star. This is because along the same LOS path the red- and 
blue-shifted masers would appear in front of and behind the star, respectively, as long as the 
coherence path lengths satisfy the requirement of maser excitation, under the condition that 
the maser gas was undergoing infall to the central star during our observations. However the 
blue-shifted masers generated behind star would be obscured by the stellar disc projected on 
the LOS, while the red-shifted masers would not. Thus only red-shifted maser emission will 
be seen closer to the center. Actually, we have confirmed that the SiO maser shell contracts 
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to the star during our observations in Sect. 3.2.2. Moreover, from Figure 3, some systemic 
masers with velocity —7 km s -1 (see Sect. 3.2.1) locate at the distance of 7 mas which can 
be viewed as an upper limit to the photospheric radius, then the extreme blue-shifted masers 
with the distance of less than 7 mas would be obscured. This is consistent with our data 
(see Figure 3) showing that all blue-shifted masers locate at the distance of larger than 7 
mas. 



From Figure 2, we notice that there appears a velocity gradient at both epochs, with the 
bluest- and reddest-shifted maser features lying closer to the center of the distribution than 
those with intermediate velocities. To verify this, we plotted in Figure 3 feature distance 
from the fitted center (marked by the red star in Figure 2) versus its LOS velocity for both 
epochs. Both epochs appear to have the same distribution with a peak near the velocity 
of —7 km s -1 and decreasing maser distance with increasing deviation of velocity from this 
peak velocity. This has also been seen in some OH maser sources (e.g. Reid et al. 1977, 
Chapman & Cohen 1986), H2O maser sources (e.g. Yates & Cohen 1994) and SiO maser 
sources (e.g. Boboltz & Marvel 2000; Wittkowski et al. 2007). A widely used simple model 
to explain this phenomenon is that of a uniformly expanding thin shell (e.g. Reid et al. 
1977; Yates & Cohen 1994; Wittkowski et al. 2007). In this model the projected distance r 
(in the fourth column of Tables 1 and 2) of a maser on the shell from the center is related 
to its LOS velocity Vlsr by the expression 



where V* is the systemic velocity of maser source, r s is the shell radius and V errp is the 
expanding velocity. Apparently, this model traces an ellipse on the r — Vlsr plot. 

The uniformly expanding thin shell model was used to characterize the expansion or 
contraction kinematics of a circular maser distribution. For the case of AH Sco, even though 
we characterize the maser distribution as an ellipse, an axial ratio of ~1.2 suggests that the 
distribution of maser features is approximately a circular structure as discussed in Section 
3.1. Thus, we can also apply the uniformly expanding thin shell model to AH Sco. Moreover, 
most maser features locate in the northwest and southeast (i.e. the direction of major axis), 
and only few maser features locate in the northeast and southwest (i.e. the direction of minor 



3.2. 



The kinematics of the SiO masers 



3.2.1. The kinematics obtained from LOS velocities 
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axis). These make the assumption of uniformly expansion/contraction in all directions of the 
uniformly expanding thin shell model to be suitable for the case of AH Sco. We performed 
a least-squares fit of the uniformly expanding thin shell model to the distribution of epochs 
A and B. The LSR stellar velocity of AH Sco, which is estimated to be about —7 and —3 
km s _1 based on the observations of OH maser (Baudry, Le Squeren & Lepine 1977) and 
H2O maser (Lepine, Pase de Barros & Gammon 1976), has not been measured particularly 
well yet. Thus the LSR stellar velocity V*, together with the shell radius r and expansion 
velocity Y exp , is treated as a free parameter in the fitting procedure. The values of V*, Y exp 
and r s were found to be — 6.8±0.5 km s _1 , 18.8±2.0 km s _1 and 9.3±0.1 mas, respectively, 
for both epochs, where the uncertainties are their standard errors. The best fitted r — Vlsr 
ellipse is also plotted in Figure 3. The systemic velocity of AH Sco of —6.8 km s _1 is 
consistent with the value of —7 km s _1 from the OH maser observations. The shell radius of 
9.3 mas determined from above model fitting is larger than the radius of maser distribution 
of 8.6 mas (see Sect. 3.1). A note is that the definition of the shell radius determined from 
the uniformly expanding thin shell model is different from that of the radius of the maser 
distribution. The shell radius reflects the scale of a 3-dimensional maser spherical shell, 
whereas the radius of maser distribution reflects the scale of maser distribution in the sky 
plane, and is the projected size on the sky plane of the 3-dimensional spherical shell. Thus 
it is not surprising that the shell radius is a bit larger than the radius of maser distribution. 

For comparison, the escape velocity calculated at the shell radius of 9.3 mas, assuming 
a typical mass of 10 M Q for supergiant and a distance to AH Sco of 2.26 kpc, is about 29 
km s -1 . Thus, the expansion/contraction velocity at the location of SiO maser shell is less 
than the corresponding escape velocity, suggesting that the maser gas is still gravitational 
bound to the star. However, the current fitting can not tell the sign of V exp term (as can be 
seen in Eq. (1)), and thus can not differentiate between expansion and contraction of the 
maser shell. The dominant expansion or contraction of maser shell can be clarified by the 
SiO maser proper motion analysis to be discussed below. 



3.2.2. Maser Proper Motions 

We can study proper motion and the kinematics of the CSE of AH Sco by tracing the 
matched features that appeared in both epochs. Because the absolute position of the phase 
center in each image is not kept during the data reduction, we must align two-epoch maps 
for studying the proper motion. The feature used for registration is the one with a velocity 
Vlsr ~ —9.3 km s _1 (labelled by "R" in Tables 1 and 2) at both epochs. And then we shift 
the coordinate frames for both epochs to align the origin (0, 0) with this feature. At an 
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assumed distance of 2.6 kpc (Baudry et al. 1977) and a maximum expansion/contraction 
velocity of 20 km s -1 (see Sect. 3.2.1), the maser proper motion should be less than 0.05 
mas in an interval of 12 days. Thus we can match these features from one epoch to another 
epoch using the criterion that the angular separation of the matched features between two 
epochs should not exceed 0.15 mas after allowing the maximum position uncertainty of (0.05, 
0.10) mas (see Tables 1 and 2) and they have similar velocity profile and flux density. As a 
result, we identified 59 commonly matched maser features between two epochs (see Tables 
1 and 2). Actually, using a reference feature located on the maser shell to align maps for 
two epochs could introduce a constant offset vector representing the motion of the reference 
feature in the individual maser proper motions. We assumed that the vector-average of the 
proper motions for all the matched features represents the motion of the aligned feature. 
In order to present a better representation of the real motions of individual features, the 
mean proper motion was subtracted from each of the determined proper motion vectors. 
The proper motions of matched maser features are shown in Figure 4. Here we adopt the 
distance to AH Sco of 2.26 kpc (see Sect. 3.2.3) for denoting the velocity values of the proper 
motions. 

From Figure 4, we can clearly see that the maser shell shows an overall contraction 
toward the central star. In order to better characterize the net contraction of the masers, 
we computed the separations between pairwise combinations of features. This technique 
has previously been applied to analyse proper motions of OH masers (Chapman, Cohen & 
Saika 1991; Bloemhof, Reid & Moran 1992), H 2 masers (Boboltz & Marvel 2007), and SiO 
masers (Boboltz, Diamond & Kemball 1997; Chen et al. 2006), and has no dependence on 
the alignment of maps. The procedure involves computing the angular separation between 
two features at the first epoch and the separation between the corresponding two features 
at the second epoch. The difference between the two values of separation is referred to as 
the pairwise separation. The procedure is repeated for all the possible pair combinations. 
However, the inclusion of all possible pair combinations often results in decreasing toward 
zero due to the bias caused by calculating pairs of closely spaced maser features. For the sake 
of clarity, and to determine representative values for the angular shifts due to the contraction, 
we have included only those pairs separated by more than 9 mas (corresponding to the 
radius of maser distribution). We obtained the mean value of these pairwise separations 
of —0.039 ± 0.002 mas, in an interval of 12 days, corresponding to a proper motion of 
— 1.186 ±0.061 mas yr _1 or a velocity of —12.7 ±0.7 km s _1 at a distance of 2.26 kpc, where 
the uncertainties are the standard errors. The negative value of proper motion implies an 
overall contraction of the maser shell. The contraction value of 1.186 ± 0.061 mas yr _1 of 
maser shell derived from the pairwise separation is significantly less than the scalar-averaged 
value of SiO proper motions of 1.96 ± 0.15 mas yr _1 (where its uncertainty is the standard 
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error). This is because some of maser proper motions do not completely point to the center 
or even a few proper motions show outflow motion as can be seen in Figure 4. 

The contraction of SiO maser shell has been reported in two Mira variables R Aqr 
(Boboltz, Diamond & Kemball 1997) and TX Cam (Diamond & Kemball 2003) and one red 
supergiant VX Sgr (Chen et al. 2006; 2007). Our observations provide an inward motion 
of SiO maser shell around another red supergiant AH Sco. In Table 3 we list these four 
sources. We also estimated the stellar optical phase of AH Sco at our observation sessions 
to be ~ 0.55 (i.e. at the optical minimum phase) based on the American Association of 
Variable Star Observers (AAVSO) data. Interestingly, the optical phase at which the SiO 
maser shell around the red supergiant AH Sco contracts is nearly the same as that seen in 
other three sources: VX Sgr (0 = 0.75 - 0.80; Chen et al. 2006), R Aqr (0 = 0.78 - 0.04; 
Boboltz, Diamond & Kemball 1997), TX Cam (0 = 0.50-0.65; Diamond & Kemball 2003). 
This infers that the contraction of the SiO maser shell would occur during an optical stellar 
phase of 0.5 — 1, which agrees with the previous conclusion reported by Chen et al. (2006) 
and the theoretical kinematical model results of Humphreys et al. (2002). Moreover, from 
Table 3 we can find that the contraction velocity of about 13 km s _1 of maser shell around 
AH Sco is the largest among the four sources. 



In order to estimate further the kinematical parameters of SiO masers and the distance 
to AH Sco, we made model-fitting to analyze spatial distribution and proper motion of SiO 
maser features as done by Gwinn, Moran & Reid (1992) and Imai et al. (2000; 2003). The 
model fitting requires to minimize the squared sum of the differences between the observed 
and model velocities, 



where, fii X and are the observed proper motions in R.A. and Dec, respectively, and Ui Z 
the observed velocity along LOS, d the distance to the maser source, a = 4.74 km s _1 mas 
_1 yr kpc -1 , (<Ji X , <Jiy, <Ji Z ) the standard deviations of the observed velocity vectors which are 
determined in the similar manner of Imai et al. (2002). In this work, we assume a spherically 
expanding flow in SiO maser region, thus the model velocity vector Vj (including V^, Vi y in 
R.A. and Dec, and V{ Z in LOS) for the ith maser feature can be expressed as 



3.2.3. 3- dimensional kinematics model for SiO masers 



x 2 = £{ 




(2) 



(3) 



I 
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Ti = Xi- x (or r ix = Xi - x , r iy = y t - y , r iz = zj, (4) 

where V (v 0x , v 0y , v 0z ) is a systemic velocity vector of the stellar system, reflecting the 
motion of the central star; V exp (i) is an expanding velocity as a function of the distance from 
the origin of the flow, r,; x (xo, yo, 0) is the position vector of the flow origin with respect 
to the position of reference maser feature (here, we assumed that the origin of the flow was 
at the star, whose positions of Xq = —7.94 mas and y = —5.83 mas have been derived from 
the least-squares fit of an ellipse to the distribution of masers in Sect. 3.1.); and (xi, yi) is 
an observed position of maser feature on the sky plane; the position of a maser feature along 
the LOS, Zi is estimated as one of the free parameters too (Imai et al. 2000). In the model 
fitting procedure, we adopted the expanding velocity as V exp {i) = Vi(rj/r ) a , where V\ is the 
expansion velocity at a unit distance of r = 10 mas, and a a power-law index indicating 
the apparent acceleration or deceleration of the flow. Moveover, we excluded those maser 
features with large positive expansion velocities in the fitting. Finally, we used 48 proper 
motion data and obtained the best solutions with their standard errors, which are given in 
Table 4. 

From Table 4, we can see that a negative expansion velocity of —14.1 ± 1.4 km s _1 
at a shell radius of 10 mas was estimated from the best-fit model, supporting the presence 
of a real contracting flow in the SiO maser region around AH Sco. This is consistent with 
the conclusion derived from the pairwise separation calculation in Sect. 3.2.2. For com- 
parison, we also estimate the corresponding velocity of —15.0 ± 1.5 km s _1 at the mean 
SiO maser radius of 9 mas according to the power-law index of the acceleration of the flow 
a obtained in our model. However, the velocity of —15.0 ± 1.5 km s" 1 seems larger than 
that of —12.7 ± 0.7 km s~ x obtained from the pairwise separation analysis, which is because 
that we excluded maser features showing outflow motion in the 3-dimensional kinematics 
model fitting procedure. And the velocity of 15.0 ± 1.5 km s~ x is slightly smaller than that 
of 18.8±2.0 km s _1 obtained from a least-squares fit of the uniformly expanding thin shell 
model to SiO maser LOS velocity structure (see Sect. 3.2.1), which is due to that some 
maser proper motions used in model fitting still deviate from the originating point of the 
inflow (i.e. the position of central star; see Figure 4), whereas the model fitting method 
involves one critical assumption that velocity vectors are in the radial direction from the 
commonly originating point of the flow. The systemic velocity along LOS of v 0z = —5.2 km 
s _1 from the 3-dimensional kinematics model is also roughly consistent with that of —6.8 km 
s _1 derived from the uniformly expanding thin shell model in Sect. 3.2.1, suggesting that 
the assumption of the origin of the flow located at the star is reasonable. More interestingly, 
we obtain a negative power-law index of the acceleration of the flow, a = —0.54 ± 0.16, 
indicating that the contracting flow was accelerating in the SiO maser region with a similar 
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form of v ~ r~ 5 of the gravitational contraction. Thus the 3-dimensional maser kinematics 
model suggests that the infall motions of SiO masers can be achieved under the gravitational 
effect of the central star. 

The velocity gradient across the SiO maser region usually used in the numerical sim- 
ulation of SiO masers (e.g. Doel et al. 1995; Humphreys et al. 2002) can be expressed 
by 

_ dlnV _ rdV 
£ ~ dim ~ Vd? ^ 
A value of e = corresponds to a constant velocity expansion, while e > 1 corresponds 
to a velocity field with large radial accelerations. A value of e — 1 was usually adopted in 
the SiO maser numerical simulations (e.g. Doel et al. 1995). Chapman & Cohen (1986) 
estimated the velocity gradient e for the OH, H 2 and SiO maser emission around VX Sgr 
and found that e ~ 0.2 in the 1612 MHz OH maser region, e ps 0.5 in the region of the H 2 
and mainline OH masers, and e ~ 1 in SiO maser region. However, the power-law index of 
the acceleration of the flow (a = —0.54 ± 0.16) derived from our best-fit kinematic model 
suggests a negative velocity gradient value of e — —0.54 in SiO maser region around AH 
Sco. This is different from those flows that apparently exhibit the accelerations in their SiO 
maser kinematics (e.g. VX Sgr, Chapman & Cohen 1986; S Ori, Wittkowski et al. 2007) 
and the positive velocity gradient value used in the maser simulations. The content of maser 
simulations is beyond this work. However, we think that such a negative velocity gradient 
e adpoted in maser simulation may be necessary for understanding the SiO maser emission, 
especially during the infall stage of the SiO maser shell. 

The 3-dimensional kinematics model fitting shows a best solution for the distance to 
AH Sco of 2.26 ± 0.19 kpc. This distance value seems reasonable. Firstly, this distance of 
AH Sco is in a good agreement with its estimated 'near' kinematic distance of about 2.0 
kpc at the systemic velocity of —6.8 km s _1 with the adopted galactic constants, R = 
8.5 kpc and O = 220 km s^ 1 . Secondly, at the distance of 2.26 kpc, the scalar-averaged 
value of 1.96 ± 0.15 mas yr _1 of proper motions of matched features (shown in Figure 4) 
would correspond to a velocity of 21.0 ± 1.6 km s _1 . This velocity is consistent with the 
expansion/contraction velocity of 18.8 ± 2.0 km s -1 in SiO maser region obtained from a 
least-squares fit of the uniformly expanding thin shell model to LOS velocity structure. Thus 
we adopt the distance to AH Sco of 2.26 kpc throughout this work. 



- 12 - 



4. Conclusions 

We summarize the main results obtained from 2-epoch (at an interval of 12 days) mon- 
itoring observations of the 43 GHz v = 1, J = 1 — SiO maser emission toward AH Sco 
performed in March 2004, corresponding to a stellar optical phase of ~ 0.55. 

(1) . Our observations revealed a persistent elliptical structure of SiO masers with the sizes 

of the major and minor axes of about 18.5 and 15.8 mas, respectively, along a position 
angle of 150°. We notice that the red-shifted SiO maser emission lies slightly closer to 
the center than the blue-shifted one. 

(2) . The LOS velocity structure of the SiO masers shows a velocity gradient at both epochs, 

with masers towrad the blue- and red-shifted ends of the spectrum lying closer to the 
center of the maser distribution than masers at intermediate velocities, which can be 
explained by the outflow or infall kinematics of SiO maser shell. By analyzing the 
uniformly expanding thin shell model to the LOS velocity of SiO masers, we estimated 
the expansion/contraction velocity of about 19 km s" 1 in SiO maser region around AH 
Sco. 

(3) . The proper motions of 59 matched features between two epochs show that the SiO 

maser shell around AH Sco was undergoing inward motion to the central star. Com- 
puting pairwise separation of these matched features, we obtained that the maser shell 
contracts toward AH Sco with a velocity of about 13 km s _1 at a distance to AH Sco 
of 2.26 kpc. The stellar optical phase of red supergiant AH Sco is very close to that of 
Mira variables R Aqr and TX Cam and red supergiant VX Sgr when the SiO maser 
shell contracts. And the contraction velocity of about 13 km s _1 of maser shell around 
AH Sco is the largest one among the known four sources showing contraction of SiO 
maser shell. 

(4) . We made a 3-dimensional kinematics model to analyze spatial distribution and proper 

motion of SiO maser features. The 3-dimensional maser kinematics model further 
suggested that the contraction of SiO maser shell around AH Sco is mainly due to the 
gravitation of the central star. And the distance to AH Sco of 2.26±0.19 kpc estimated 
from this kinematics model fitting is consistent with the kinematic distance of 2.0 kpc 
at the systemic velocity of AH Sco of —7 km s -1 . 
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Table 1. 43 GHz SiO maser features around AH Sco observed by VLBA on March 8, 2004. 
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° 'x 
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/) 
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11 


,7 
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8 


1 
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6.96 


-1 


.356 
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,7 
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19 


,7 
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5 


2 
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8.76 


-1 
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.007 
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4, 
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5 
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2 
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-1 
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16 


.6 
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3 


6 


6 
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8 
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. (1(12 
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9 


.4 
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1 1 
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2 


.17 


9.47 


-0 


.330 
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9 


.8 
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9 


13 


12 
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1 
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-1 
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4 


.5 
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9 


14 


13 
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2 
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-0 
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.9 
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8 


16 
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x 


20 
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I 
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1 


95 


9.39 


-15 


.667 





.006 


11.237 





.018 
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5 


21 
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41 
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g 
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Table 1 — Continued 



ID 
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o 'x 
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Match ID 




(km s" 


-1) 


(mas) 
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(mas 


) 




(Jy) 
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(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 




(8) 


(9) 


(10) 


(11) 


66 


4.10 


3.04 


8.66 


0.508 


0.008 


4.077 





.013 


2.8 


29.4 


70 


67 


4.70 


1.52 


7.40 


-12.565 


0.009 


11.650 





.020 


5.8 


30.0 


71 


68 


4.83 


1.30 


6.67 


-9.377 


0.012 


-0.696 





.033 


2.3 


10.2 




69 


5.20 


1.95 


7.22 


-3.566 


0.014 


0.108 





.027 


2.4 


14.6 




70 


5.43 


1.30 


6.05 


-10.355 


0.003 


11.399 





.012 


5.5 


19.7 


75 


71* 


5.59 


0.65 


6.82 


-4.651 


0.008 
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.023 


1.2 


3.2 


76 


72 


5.97 


1.52 


5.75 


-13.283 


0.009 


3.593 





.016 


2.5 


15.5 


74 


73 


7.00 


1.30 


6.21 


-7.430 


0.011 


12.023 





.025 


1.9 


8.1 


79 


74 


7.33 


3.04 


6.39 


-9.017 


0.006 


-0.477 





.013 


4.0 


41.0 


78 


75 


7.42 


1.74 


6.51 


-8.278 


0.008 


12.338 





.022 


2.9 


16.6 




76 


7.73 


2.17 


5.08 


-12.707 


0.005 


3.946 





.008 


4.0 


27.5 


82 


77* 


7.99 


0.65 


7.79 


-0.359 


0.017 


4.208 





.034 


1.3 


3.5 




78* 


8.01 


1.09 


6.16 


-8.120 


0.008 


11.995 





.023 


2.3 


7.8 


83 


79* 


8.35 


1.09 


3.98 


-4.070 


0.014 


6.607 





.026 


1.5 


6.5 




80 


8.92 


1.95 


5.85 


-8.965 


0.006 


0.059 





.023 


4.9 


31.9 


84 


81 


9.17 


1.30 


5.85 


-8.289 


0.007 


11.678 





.013 


3.9 


17.9 


85 


82 


12.62 


0.87 


2.50 


-5.479 


0.009 


5.706 





.017 


2.2 


7.4 


86 



Note. — column (l): ID number; columns (2): Vlsr a t the peak of velocity profile of feature; column (3): the velocity 
range across the feature An; column (4): distance of mascr feature, r, from the fitted position of central star; columns (5) 
and (7) : the intensity weighted ccntroid of each feature (x, y); columns (6 J and (8 J : the corresponding uncertainties (<r x , <jy ); 
column (9): the peak flux density of each feature P; column (10): the integrated flux density of all spots in the feature S; 
and column (11): the ID numbers of matched features at another epoch. 

* Feature which can not be well represented by a Gaussian curve. 

R Reference feature. 
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Table 2. 43 GHz SiO maser features around AH Sco observed by VLBA on March 20, 2004. 



ID 
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(km s" 
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(mas 
(7) 
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(8) 


P 
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Epoch 2 
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1 
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1 


09 
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-2 
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003 


-0 
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.014 


5 


.8 


16. 





1 


2 


-15.14 


2 


.39 


8.33 


-1 


.650 
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.009 


14 


,1 


90 





3 


3* 
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.87 


7.97 


-1 
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.021 


2 


.5 


3. 
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4 
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.87 
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-1 
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1 
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.011 


1 


.3 


3. 
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5 
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8.90 


-1 
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010 


-0 
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.011 


1 


,2 


3. 
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6 
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2 


.39 


8.39 


-1 
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.009 


21 


,7 


121. 


4 


5 


7 


-13. 74 
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9.04 


-0 
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.091 


2 


,1 


5. 
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8 


-13.31 





.65 


8.37 


-3 


.046 





. 005 


-0 
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.017 


2 


.0 


4. 


4 


7 


9 
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4 


.13 


8.85 


-0 


.813 





006 
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.015 


18 


9 


144. 
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9 


1(1 


-12.17 
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.82 
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Table 2 — Continued 



ID VlSR a « r x V "v p S Match ID 

(km s (mas) (mas) (mas) (Jy) (Jy km s *) Epoch 2 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) 



66 


3 


50 


1 


52 


6 


.85 


-10 


.396 





Oil 


12 
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5 


.1 


23 


.6 


62 


67 


4. 
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.283 
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.026 


2 


.5 


14 


.7 


66 


71 


4. 


78 


1 


52 


7 


.37 


-12 


.552 





013 


11. 


604 





.020 


8 


.0 


35 


.3 


67 


72 


5 


04 


1 


74 


6 


.40 


-1 


.921 





007 


8 


073 





.026 


2 


.5 


13 


.7 




73* 


5. 


44 





.87 


6 


.93 


-3 


.062 





Oil 





916 





.038 


1 


.2 


3 


.6 




74 


5. 


56 


1 


30 


5 


.79 


-13 


.234 





012 


3 


534 
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.87 


6 


.84 


-4 


.629 





011 


-0 


132 





.038 


1 


.5 


5 


.7 


71 


77 


6. 


06 





.87 


7 


13 


-3 


.665 





008 





137 





.026 


1 


.7 


5 


.4 




78 


7. 


05 


1 


.74 


6 


.46 


-8 


.987 





013 


-0 


517 





.027 


5 


.3 


25 
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84 
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Note. — The representations of columns (l)-(ll) are the same as in Table 1. 
* Feature which can not be well represented by a Gaussian curve. 
^Reference feature. 
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Table 3. Sources with the detected contraction of SiO maser shell. 



Source 


stellar phase 


Contraction velocity 


Reference 







km s _1 




R Aqr 


0.78 - 1.04 


4.2±0.9 


Boboltz et al. (1997) 


TX Cam 


0.50-0.65 


5-10 


Diamond & Kemball (2003) 


VX Sgr 


0.75 - 0.80 


4.1±0.6 


Chen et al. (2006) 


AH Sco 


~ 0.55 


12.7±0.7 


this work 



Table 4. Best-fit model for the SiO maser kinematics in AH Sco. 



VQx 


V y 


V0z 


Vi 


a d 


Reduced % 2 


(km s _1 ) 


(km s _1 ) 


(km s _1 ) 


(km s _1 ) 


(kpc) 




3.8 ±0.3 


-0.9 ±0.4 


-5.2 ±0.4 


-14.1 ± 1.4 


-0.54 ±0.16 2.26 ±0.19 


3.58 
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Fig. 1. — Top: Comparison of total power (solid line) to cross power (open circle) of 43 
GHz v=l, J— 1-0 SiO maser emission toward AH Sco obtained on (a) March 8, 2004 and (b) 
March 20, 2004. Bottom: The corresponding fraction power (cross/total) detected by the 
high-resolution VLBA observations. 
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Fig. 2. — VLBI images of 43 GHz v=l, ,7=1-0 SiO maser emission toward AH Sco obtained 
on (a) March 8, 2004 and (b) March 20, 2004. Each maser feature is represented by a filled 
circle whose area is proportional to the logarithm of the flux density, and the color indicates 
its Doppler velocity with respect to the local standard of rest. Its stellar velocity is about 
—7 km s^ 1 . Errors in the positions of the features are smaller than the data points. The 
ellipse indicates the least-squares fit to the maser distribution for each epoch. The fitted 
center of ellipse model is marked by the red star. 
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Fig. 3. — Distance of maser features from the fitted position of central star (in Fig. 2) versus 
their LOS velocity for epochs A and B. Maser features of the different epochs are denoted by 
different color symbols whose area is proportional to the logarithm of the flux density. The 
plot suggests that the higher-velocity maser features lie closer to the central star, which can 
be well explained by the uniformly expanding thin shell model. Indicated by the downward 
arrow is the systemic velocity of AH Sco of —6.8 km s _1 , obtained from the best-fitting 
thin-shell model (shown by the curve). 
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Fig. 4. — Distribution of proper motion velocity vectors of the matched maser features at 
an assumed distance of 2.26 kpc. The length of the vector is proportional to the velocity. 
The mean proper motion vector has been subtracted from each of the determined proper 
motion vectors. The color and size of symbols are the same as that shown in Fig. 2. Red 
star represents the fitted center of ellipse model to maser distribution (see Sect. 3.1). 



